Ground states of a frustrated quantum spin chain with long-range interactions 
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The ground state of a spin-1/2 Heisenberg chain with both frustration and long-range interac- 
tions is studied using Lanczos exact diagonalization. The evolution of the well known dimerization 
transition of the system with short-range frustrated interactions (the J1-J2 chain) is investigated in 
the presence of additional unfrustrated interactions decaying with distance as l/r a . It is shown that 
the continuous (infinite-order) dimerization transition develops into a first-order transition between 
a long-range ordered antiferromagnetic state and a state with coexisting dimerization and critical 
spin correlations at wave-number k — tt/2. The relevance of the model to real systems is discussed. 

PACS numbers: 75.10.Jm, 75.10.Nr, 75.40.Mg, 75.40.Cx 



One-dimensional spin systems have played an impor- 
tant role in quantum many-body physics since the early 
days of quantum mechanics [H 0] . Several different types 
of ordered and disordered ground states can be realized, 
depending on the individual spin magnitude S and the 
form of the spin-spin interactions [3H6[. For S = 1/2, 
the prototypical Heisenberg chain with antiferromagnetic 
nearest-neighbor interactions (coupling constant J\ > 0) 
has a quasi-ordered (critical) ground state, with spin cor- 
relations decaying with distance r as ydrr(r)/r |fjj. In- 
cluding a next-nearest-neighbor coupling J2 > (the Ji- 
J2 chain Q) leads to a quantum phase transition into a 
doubly-degenerate dimerized state (a valence-bond-solid; 
VBS) at coupling ratio g = J 2 /J\ ~ 0.2411. In the ef- 
fective field theory for the S = 1/2 chain @, the VBS 
transition is related to a sign change of a marginal oper- 
ator. It has been investigated in great detail numerically, 
using, e.g., exact diagonalization 0, Q and the density- 
matrix renormalization-group (DMRG) method [9]. 

While long-range spin ordering is rigorously ruled out 
in one-dimensional systems with finite-range rotationally 
invariant interactions, long-range interactions make mag- 
netic order possible at zero temperature. The transition 
between a long-range ordered antiferromagnet (AFM) 
and the quasi-long-range ordered (QLRO) ground state 
was recently investigated in a Heisenberg chain with in- 
teractions of the form J r oc (— l) r ~ 1 /r a [10j. Here the 
signs correspond to no magnetic frustration, thus favor- 
ing AFM ordering. For a < er c , the ground state pos- 
sesses true AFM long-range order, while for a > a c the 
system is in a QLRO phase, with the same critical form 
of the spin correlations as in the standard Heisenberg 
chain. The critical value ev c depends on details of the 
couplings (e.g., on J\ when all other J r are fixed) and 
the exponents are continuously varying. 

Another example of long-range interactions is the cel- 
ebrated Haldane-Shastry chain ll|, with frustrated in- 
teractions J r — 1/r 2 . This system has a critical ground 
state similar to that of the standard Heisenberg chain, 
but the marginal operator vanishes [12| and it is, thus, a 
system right at the dimerization transition. 



A natural question arising from previous work is how 
the combined effects of frustration and long-range inter- 
actions could lead to other phases and quantum phase 
transitions. In particular, is it possible to realize a di- 
rect transition between the AFM state and a VBS? In 
this Letter the evolution of the standard dimerization 
transition into an AFM-VBS transition is explored by 
considering a frustrated J1-J2 chain with additional non- 
frustrated long-range interactions. The hamiltonian for a 
finite periodic chain with N spins S = 1/2 is given by 



N/2 N 

H = J r Si ■ Sj_|_ r , 

r=l i=l 

where the couplings are given by 



■h = g, Jr^2 = 



(1) 




(2) 



Here the normalization is chosen such that the sum of all 
non- frustrated interactions \J r \ equals 1 [l~3| (and J\ is 
also given by the J r ^2 expression). 

The model is here studied using Lanczos exact diag- 
onalization. A semi-quantitative phase diagram based 
on these calculations in the plane {g,a~ l ) is shown in 
Fig. [2 The J1-J2 chain corresponds to the horizontal axis 
{a- 1 = 0). The QLRO phase is here denoted QLRO(tt), 
with 7r indicating the wave-number of the dominant spin 
correlations. The phase boundaries are approximate, re- 
sulting primarily from studies of level crossings, as will 
be discussed below. The main focus of this initial study 
of the model will be on the evolution of the QLRO(7r)- 
VBS (dimerization) transition with decreasing a. It will 
be shown that this continuous transition persists until 
a r; 2, while for smaller a it evolves into a first-order 
transition (of the avoided level-crossing type) between 
the AFM state and a state with coexisting VBS order 
and critical spin correlations at wave-number k = tt/2, 
denoted in the phase diagram as VBS+QLRO(7r/2). 

The coexistence state is not purely of theoretical in- 
terest. Recent ab-initio calculations for metallic chains 
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FIG. 1: (Color online) Approximate phase diagram as a func- 
tion of the frustration strength g and the inverse of the long- 
range exponent a. The dashed curves indicate continuous 
transitions, whereas the thick solid curve represents a first- 
order transition. The curve for g < corresponds to the 
interaction studied in At a -1 = the dominant spin 

correlations in the VBS state change from k = tt to n/2 at 
g re 0.52 This transition (or cross-over) evolves to the 
point where all phase boundaries come together. 

show unfrustrated spin couplings decaying as re 1/f* 2 , 
with Ji in some cases frustrating (e.g., Mn) |l4[. In the 
quasi-classical (large-S") limit, spiral states with contin- 
uously varying periodicity can arise in such a systems. 
The present study suggests a more exotic scenario in the 
extreme quantum limit of S — 1/2 (and perhaps also for 
other small S). 

'Solving the model (1) numerically poses significant 
technical challenges. Efficient quantum Monte Carlo 
techniques can be applied to systems with long-range 
interactions [HI [H|, but with the frustrating J2 term 
this is no longer possible due to the sign problem [lij . 
The DMRG method [H, on the other hand, can 
handle frustration but not easily long-range interactions. 
Here periodic chains up to size N — 32 are solved using 
Lanczos exact diagonalization (in the standard way, ex- 
ploiting lattice symmetries and spin-inversion for block- 
diagonalization in the magnetization m z — sector). 
This is sufficient for roughly extracting the phase bound- 
aries using level crossing methods (which in the case of 
the dimerization transition is a well established technique 
0, extended here using different levels to detect other 
transitions) . 

The QLRO(-7r)-VBS transition in the J1-J2 chain is of 
infinite order, i.e., the singlet-triplet gap of the VBS is 
exponentially small for .9 — > g c [5J- ^ is therefore diffi- 
cult to locate the transition based on the order parameter 
for small N. However, g c can be determined accurately 
from excited states. The lowest excitation of a chain with 
even N is a triplet for g < g c and a singlet for g > g c . 
The crossing point of these levels is a rapidly converg- 
ing finite- N definition of g c 0, H|- The same physics 
can be expected also in the presence of the long-range 
interaction, if a is sufficiently large. This is shown for 
a 16-spin chain at a = 4 in the upper panel of Fig. [2] 
Singlets with momenta k = and k = n (out of which 
symmetry-broken dimerized states can be formed) should 
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FIG. 2: (Color online) Low-energy levels of a 16-spin system 
at q = 4 (upper panel) and a = 1 (lower panel). The spin S 
and the momentum k of the states are indicated in the upper 
panel. The inset in the lower panel shows the avoided level 
crossing of the two k — singlets in greater detail. 

be degenerate in the VBS phase. For finite ./V this de- 
generacy is not exact (except in the J1-J2 chain at the 
special point g = 1/2), but a region of very near degen- 
eracy for g > 1/2 can be seen in the figure. The region of 
approximate degeneracy, which is not easy to demarcate 
precisely, expands very slowly toward smaller g with in- 
creasing N. In contrast, the singlet-triplet crossing point 
is well defined and converges rapidly. Extrapolating the 
crossing point to N — 00 for different a, as illustrated in 
Fig- El can reliably give the QLRO(7r)-VBS phase bound- 
ary g c (p) for a > 2. 

Upon decreasing a below re 2, the broad maximum 
in the ground state energy versus g becomes increasingly 
sharp. As seen in the lower panel of Fig. [3J at a = 1 it has 
developed into a sharp tip due to an avoided level crossing 
with the second singlet at k = 0. The real singlet-triplet 
crossing has moved to the same region. An avoided level 
crossing leading to a discontinuity in the derivative of 
the ground state energy with respect to g for N — > 00 
is the hall-mark of a first-order transition. The nature 
of the phases at this transition will be discussed below. 
First, let us investigate how the transition evolves from 
continuous to first-order. 

Fig. [3] shows the size dependence of the level crossing 
point Across and the location (? pca k of the maximum in the 
ground state energy. In the Ji- J2 chain the size correction 
to the crossing point is cx 1/JV 2 , which also can be seen 
for large a. For smaller a, the corrections instead seem 
to be oc 1/iV, but a cross-over to 1/JV 2 for large N seems 
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FIG. 3: (Color online) Dependence on the inverse chain length 
of the singlet-triplet crossing point Across and the location 
Speak of the ground state energy maximum for different long- 
range interaction exponents a. The two lines show extrapo- 
lations of the g = 2.0 numerical data to N — oo. 

likely as long as the transition remains continuous. The 
peak location moves in the opposite direction. For some 
a and N — > oo, Across and <7 pea k should coincide. The 
results indicate that both g croS s and g pca k have dominant 
1/N corrections at this point. Fitted lines are shown in 
Fig. |3] at a = 2, were there is still a small gap between 
the two extrapolated values. For a — 1.7, where the 
transition is first-order, they should coincide (and then 
the asymptotic size correction should be exponential). 

To verify an avoided level crossing with a discontinuous 
energy derivative for a < 1.8, the second derivative of 
the ground state energy at its maximum is graphed on 
a lin-log scale in Fig. [4] It grows exponentially with N 
for a — 1.5, showing that the slope of the energy curve 
indeed changes discontinuously for an infinite chain. In 
contrast, at a = 3 the second derivative decreases for 
large N. For a — 2 convergence to a finite value also 
seems plausible, whereas a = 1.7 and 1.8 appear to be 
close to a separatrix (where the form of the divergence 
is consistent with a power law) between the two different 
behaviors. 

This analysis suggests that the continuous dimeriza- 
tion transition changes smoothly into a first-order tran- 
sition at (g m rj 0.41, a m k, 1.8). The singlet-triplet cross- 
ing moves toward the ground state energy maximum and 
coincides with it at the multi-critical point (g m , Oi m ), be- 
yond which it develops into a first-order singularity. Note 
that the rounded energy maximum in the VBS phase for 
large a has no special significance (except at a -1 = 
where it corresponds to the exact singlet-product ground 
state) . It is only when it develops into the sharp avoided 
level crossing that it is associated with a phase transition. 

To discuss the states involved in the first-order transi- 
tion, consider the spin and bond correlation functions; 

C(r) = (S,.S l+r ), (3) 
D(r) = ((S. i -S 4+1 )(S l+r .S l+r+1 )>. (4) 

In Fig. [5] these are graphed for two g values, at either 
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FIG. 4: (Color online) Size dependence of the second deriva- 
tive of the ground state energy with respect to the frustration 
parameter g at the point g poa k where the ground state energy 
takes its maximum value. 

side of the transition for a — 1. At g < g c the dominant 
spin correlation C(k) in Fourier space is at wave- number 
k = tt, and finite-size scaling shows that the sublatticc 
magnetization remains non-zero for N — > oo. There is 
no structure in D(r), i.e., there is no VBS order. This is 
thus an AFM phase; the continuation of the AFM state 
studied in [lOj, as indicated in Fig. [TJ For g > g c there 
is VBS order. Interestingly, in this phase there are also 
strong spin correlations at k = tt/2, which can be seen 
clearly as a real space period-four oscillation in Fig. [5] 
Finite-size scaling indicates that there is no long-range 
spin order, but the correlations appear to decay as 1/V 7 
with 7 ~ 1; thus this state is denoted as QLRO(7r/2). 

Examining the correlations as a function of g, discon- 
tinuities (increasingly sharp jumps with increasing N) 
develop for a < 1.5. This should persist until the multi- 
critical point at a m ss 1.8, but larger systems are needed 
to observe the discontinuity very close to this point. 

The VBS+QLRO(7r/2) state should have gapless spin 
excitations. The lowest triplet has k — tt/2. It is, how- 
ever, difficult to demonstrate the gaplessness based on 
data for small systems, because the size-dependence of 
the gaps (and other quantities) for N — An exhibit even- 
odd oscillations in n. In the VBS phase the lowest triplet 
is at k = tt, even when the spin correlations (exponen- 
tially decaying) are peaked at k = tt/2. The level cross- 
ing between the lowest k = tt and k — tt/2 triplets can 
be used to extract the boundary between the VBS and 
VBS+QLRO(7r/2) phases. The size dependence of the 
crossing point is not smooth, however, and cannot be 
extrapolated very reliably. The boundary between dom- 
inant k — 7r and k = -k/2 spin correlations in the VBS 
phase has also not been extracted accurately. 

Let us briefly return to Fig. [2] for another interesting 
feature of the level spectrum: The lowest singlet excita- 
tion for small g has momentum k = n for a = 4 but 
k = for a = 1. The switching of the order of these 
levels as a function of a for g < g c is associated with the 
QLRO(7r)-AFM transition. The level crossings can be 
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FIG. 5: (Color online) Spin (upper panel) and dimer (lower 
panel) correlations in a 32-spin chain at a — 1. At g = 0.25 
and 0.45 the system is in the AFM and VBS-QLRO(tt/2) 
phase, respectively. A first-order transition between these 
states occurs at g » 0.39. 

used to extract this phase boundary very accurately up 
to g » 0.25 (while for higher g the TV — > oo extrapolations 
become difficult); a more detailed discussion of this issue 
is given as a footnote As indicated in Fig. [TJ a c de- 
pends only weakly on g. The results are consistent with 
the location quoted above for the multi-critical point. 

In summary, the combination of short-range frustra- 
tion and long-range unfrustrated interactions in one di- 
mension has been shown to lead to a first-order tran- 
sition between a Neel state and a VBS with coexist- 
ing critical k = n/2 spin correlations. It should be 
noted that the system sizes studied here are small, and it 
cannot be excluded that the spin correlations could be- 
come incommensurate, as they do in the J1-J2 chain for 
J2/J1 > 1- 20J Hopefully, field-theories that very success- 
fully describe the standard dimerization transition Q, 
and recently also the transition between the critical spin 
state and the Neel state [h| , could be generalized to the 
coexistence state as well. 

Recent calculations [14J for metallic chains have shown 
that interactions of the type used here are realistic, but 
in these systems S > 1/2. Although one cannot de- 
scribe these systems completely using a spin-only model, 
it would still be interesting to repeat the calculations 
reported here for larger S. This is much more chal- 
lenging, however, because of the rapidly growing size of 
the Hilbert space with S. Although the DMRG method 
[l7L [lU is not ideally suited for systems with long-range 
interactions, it may still be possible to use it to study 
lattice sizes beyond the limits of Lanczos calculations. 
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